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I. INTRODUCTION 

The general solution of a quantum mechanical Hamil- 
tonian problem containing a time-independent potential 
gives rise to both bound and unbound eigenstates. Usu- 
ally the Hamiltonian of the system has a finite number 
of bound eigenstates while the unbound ones form a con- 
tinuum. For the description of nuclei in the stability 
valley usually only the bound eigenstates are considered. 
However, the development of radioactive nuclear beam 
facilities has allowed the study of nuclei far from the line 
of stability, bringing to the fore new nuclear structure 
problems. One of the main topics in recent years has 
been the study of halo nuclei 0,0]. These are weakly- 
bound, spatially extended systems, typically comprising 
a core and one or two valence nucleons. 

A particularly interesting example of exotic systems is 
that of Borromean nuclei, i.e., three-body composite sys- 
tems with no binary bound states. These nuclei have de- 
served special attention because their loosely-bound na- 
ture reflects a delicate interplay between two- and three- 
body forces, thus constituting a challenge to existing the- 
ories, and a motivation for the development of new ones. 
Even today, the detailed structure of the continuum spec- 
trum of these systems is not fully understood, partially 
due to the ambiguities associated with the underlying 
forces between the constituents. Due to their low binding 
energy, halo nuclei are easily broken up in the nuclear and 
Coulomb field of the target. Therefore few-body reaction 
theories, developed to extract reliable information from 
experimental data of reactions involving loosely bound 
systems, have to include, as an essential ingredient, a re- 
alistic description of the continuum part of the spectrum. 

From the theoretical point of view, the treatment of 
reactions involving loosely bound systems deals with 
the complication that breakup states are not square- 
normalizable. A convenient method to circumvent this 
problem is to replace the states in the continuum by a 
finite set of normalized states, thus providing a discrete 
basis that, hopefully, can be truncated to a small num- 
ber of states and give a good description of the contin- 



uum. Several prescriptions to construct a discrete basis 
have been proposed. For two-body composite systems, 
where true continuum states are easily calculated, one 
can use a discretization procedure in which the contin- 
uum spectrum is truncated at a maximum excitation en- 
ergy and divided into energy intervals. For each interval, 
or bin, a normalizable state is constructed by superpo- 
sition of scattering states within that bin interval. The 
method, known as the Continuum Discretized Coupled 
Channels (CDCC) method 0,0], has been very useful in 
the description of elastic and breakup observables in re- 
actions involving weakly bound two-body projectiles. An 
alternative method to obtain a discrete representation of 
the continuum spectrum is to diagonalize the two-body 
Hamiltonian in a complete set of L 2 functions such as 
Gaussians 0, 0] or Laguerre functions 0, 0]. This 
method has the appealing feature of being readily appli- 
cable to three-body systems, in which case the Hamilto- 
nian is diagonalized in a complete set of square-integrable 
functions for the 3-body Hilbert space. Several applica- 
tions of this method can be found in the literature, for 
both structure ^3] an d reaction problems [JjJ- In the 
latter case, the method constitutes a natural extension 
of the CDCC formalism for reactions with three-body 
projectiles. 

When the ground state wave function is already known, 
a useful procedure to obtain a discrete representation 
for scattering states consists of performing a Local Scale 
Transformation (LST) 0] that transforms the ground 
state wave function of the system into the ground state 
of a Harmonic Oscillator (HO) 0G3- Once the LST 
is obtained, the HO basis can be transformed by the in- 
verse LST to a discrete basis in the physical space. The 
functions in the Transformed Harmonic Oscillator (THO) 
basis are not eigenfunctions of the Hamiltonian (except 
for the ground state) but the Hamiltonian can be diago- 
nalized in an appropriate truncated basis to produce ap- 
proximate eigenvalues and eigenfunctions. This method 
has been shown to be useful for describing the two-body 
continuum in both structure 0, Q, 0] and scattering 
flfiLri"?! problems. In particular, it was shown that global 
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structure functions, related to the coupling to the con- 
tinuum, such as strength functions, are very accurately 
described using a relatively small THO basis. 



In this work we extend the THO formalism presented 
in filj. H3 | to treat the three-body continuum. In par- 
ticular we apply the method to the Borromean nucleus 
6 He. This is a very weakly bound system with a well 
developed structure, of an a cluster and two valence 
neutrons. Most of our knowledge of this nucleus comes 
from the analysis of reactions where secondary beams of 
6 He collide with stable nuclei. These experiments have 



providing a rich variety of data which can be used to 
benchmark reaction and structure models. 



Theoretically, 6 He has been the object of many stud- 
ies, using three-body methods Ejl l2fl. I2H I2H l2fll 
cluster-orbital shell models |3ll l32j . no-core microscopic 
shell models j^H and microscopic cluster mod els, and for 
various effective nucleon-nucleon interactions |2J,lJ£|]. In 
the present work, the three-body equations are solved ac- 
cording to the formalism described in Ref. j2^|. In this 
reference, several three-body methods are compared and 
applied to the ground state properties of the Borromean 
nuclei 6 He and 11 Li. One of these approaches, discussed 
in [25] , consists of a solution of the Faddeev equations in 
configuration space. In the Faddeev formalism, the total 
wave function for a three-body system is expressed as a 
superposition of three terms, one for each Jacobi configu- 
ration. This wave function, which contains both scatter- 
ing and rearrangement channels, is obtained by solving 
a set of coupled integro-differential equations. A practi- 
cal way to solve the problem is to express the three-body 
wave function in terms of hyperspherical coordinates. For 
convenience, scattering states are expanded in terms of 
a discrete basis. In Ref. [25] a family of Gauss-Laguerre 
functions was used for this purpose. In the present work, 
we choose a THO basis obtained from a 3-body calcula- 
tion of the 6 He nucleus. 



The method we present can also be used as an input 
for reaction theories. However, in the present paper we 
present only results related to the structure of 6 He. In 
particular, we investigate strength functions for several 
operators that couple the 6 He ground state to the con- 
tinuum. Results related to scattering processes induced 
by 6 He will be presented in a forthcoming paper. 



This paper is structured as follows. In Section II the 
THO formalism for a three-body system is developed. In 
Section III the formalism is then applied to the case of 
the structure of 6 He. Finally, Section IV summarizes and 
draws conclusions from this work. 




Figure 1: The Jacobi T-coordinate system used to describe 
the 6 He system. 



II. FORMALISM OF THE TRANSFORMED 
HARMONIC OSCILLATOR (THO) STATES FOR 
A THREE-BODY CONTINUUM SYSTEM 

In this Section we develop the formalism of trans- 
formed harmonic oscillator states to be applied to three- 
body continuum systems. To a large extent, this for- 
malism is a generalization of the method presented in 
For such a system we will use the hyperspheri- 
cal coordinates that are obtained from the Jacobi coordi- 
nates (see, for instance, [36]). Consider three particles of 
masses {mi, to 2 , m 3 } and whose locations are defined by 
{ri, r 2 ,r 3 }. Three sets of Jacobi coordinates {sci,^}, 
with i = 1, 2, 3, can be defined where {x{\ is the relative 
coordinate for two particles and {y^} is the coordinate 
of the third particle relative to the center of mass of the 
two particles used to define {xi}: 



X i 



Vi 




(1) 



Here \i Xi is the reduced mass of the {j, k} system and fi Vi 
is the reduced mass of the (j + k) and particle i system, 



rrijmk 
irij + m k ' 
mjirrij + m k ) 
mi + m.j + mk 



(2) 



and m is an arbitrary normalization mass that we will 
take as the nucleon mass. The indices i, j, k run through 
(1,2,3) in cyclic order. 

From these Jacobi coordinates, one can introduce the 
hyperspherical coordinates defined by the hyperradius p 
and the hyperangle a>i, 



P 

tan a; 



vl 



vl 



y! 



(3) 



with Xi = psinai and tji — pcosai. Note that the hy- 
perangle depends on the Jacobi system selected but not 
on the hyperradius. 

For the three-body system of a core, considered inert 
and spin-less, plus two valence particles, the total wave 
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function is assumed to be a product of the core intrin- 
sic wave function C (£ C ) and the two- valence-particle plus 
core relative motion wave function ^jm(^-, 2). ^jm(1 ; 2) 
depends on the relative coordinates and the spins of the 
two particles and is the solution of the Schrodinger equa- 
tion 



(T + V-E)V(1,2) = 



with 



V = Vr 



cl 



V, 



c2 



12- 



(4) 



(5) 



V C i is the interaction of the core with particle i and V12 
the interaction between particles 1 and 2. 

From this point, we work within the Jacobi T- 
coordinate system, shown in Fig. ^ an d hence all sub- 
sequent formulae refer to this coordinate set. Expanding 
using the hyperspherical harmonics (HH) re- 
ferred to these coordinates 



JM 



(1,2) =$^00 



JM 



where p is the hyperradius, = {K,l x ,l y ,l, S} labels 
each channel, Rp{p) is the hyperradial wave function, 
T Kl v (Q) is a hyperspherical harmonic in the angles f2 = 
(a, x,y), and X$(l,2) is the total spin wave function of 
particles 1 and 2. J = I + S is the total angular mo- 
mentum of the system, since the core particle is assumed 
spinless. The hyperspherical harmonics are 



T 



Khrii 



(O) = ¥* l «(a)[Y lx (x)®Y ly (y)} 



hn / 



Wp{a) = N * B (sin a)' 11 (cos a)'" 
x P J '-+ 1 / 2 ^+ 1 / 2 (cos2a), 



(7) 



(8) 



where I x and ly are the orbital angular momenta associ- 
ated with the Jacobi coordinates x and y, respectively, 
K is the hypermomentum and P£' b is a Jacobi polyno- 
mial with n = (K — l x — l y )/2. Usually the hyperradial 
part is written as Rp(p) = p~ 5 / 2 Up(p). Then, for bound 
states U(p) — > for p — > 00 with asymptotic behaviour 
U{p) -> cxp (-hp) with h 2 k 2 /2m = -E. 



A. The THO method 

The problem is now to solve the Eq. I@J and find 
its eigenvalues and eigenvectors. For that purpose, one 
needs a basis in which the Hamiltonian can be diagonal- 
ized. As mentioned above, one can expand the channel 
wave function ^pjm(p, Ft) in the HH basis as 

<W(p,n) = Rp{p)Y.(lrniSa\JM)r l * l ? mi {Si) X%. 



mi<j 



(9) 

Rp(p) is the hyperradial wave function and here we will 
use the THO method to obtain it. The basic idea of the 



THO method is to perform a local scale transformation 
(LST) of the ground state wave function of the system 
under study into the harmonic oscillator ground state 
wave function. We now obtain such a LST for a three- 
body system. The ground state of the system ^b(p, 0) 
is written as a linear combination of the basis functions, 
Eq. @, 

* B (p,n) = Y,RB^Y.^ Sa \ jBV ) T K l L^) x l- 

(10) 

Consequently, we will obtain a LST for each channel /3 
included in the ground state. This is done by transform- 
ing the hyperradial wave function for each channel in- 
cluded in the bound ground state function, Rspip), into 
the corresponding (same K) harmonic oscillator ground 
state wave function for that channel Rqr (s). In the six- 
dimensional case, these functions are 



R?k{s) = N iK s«L? + -\s*)e* V {-s*/2) (11) 



K T K+2/2 



(6) where NiK is the normalization constant 



N, 



K 



I 2r(i + l) 
T(i + K + 3) ' 



(12) 



Here Lf +2 (x) is a generalized Laguerre polynomial with 
i the order of the polynomial and which describes the hy- 
perradial excitation. The LST for a given channel /3 in- 
cluded in the bound ground state is defined by the equa- 
tion 

/ dp P b \R B (M\ 2 = / da' s' 5 |J$? ( S ')| 2 - (13) 
Jo Jo 

Once the LST for each channel pp{s) is obtained, the 
THO basis is defined by applying the inverse transfor- 
mation sp{p) to the HO wave functions generated from 
the corresponding ground state wave function, 

rW°(p) - §^RB ( P )L?+ 2 { Sf} { P n (14) 

Usually, instead of RfJ IO (p), the following hyperradial 
wave function is introduced 



u? ho (p)=p 5/2 R™°(p) 

that fulfils the orthonormality relationship 



dpU?»°(p)U™°(p) = 6 i . 



(15) 



(16) 



Thus, we obtain the THO basis 

*£52r(P,0) = p- 5/2 U^ ho (p) (17) 
x ^(l m Sv\JM)T l & mi M x s-(l*) 



It is easy to show that these states form a complete or- 
thonormal set. 
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Since the THO basis is infinite, one diagonalizes the 
Hamiltonian in a finite truncation, obtaining eigenstates 



state components l x = l y . 



*nJM{p,V) = Ye,, 



if} 



if3 



THO 



nJM^iflJM 



(19) 



Convergence of the results with the basis truncation must 
then be checked. 

The information available on the different channels in- 
cluded in the ground state wave function allows one to 
construct the corresponding LST's directly, Eq. Ill -ill . For 
channels not included in the ground state, as a general 
rule, information from one of the known (ground state) 
channels with the closest quantum labels to the channel 
of interest is used to construct the LST. One important 
point concerns the label K which governs the p K behav- 
ior of the hyperradial wave function close to the origin. 
In order to keep this behavior correct we always select 
a channel from the ground state wave function with the 
same K as the channel under study. If this is not possi- 
ble, a channel with K — 1 is used and the corresponding 
hyperradial wave function is then multiplied by p. 

So, for example, the hyperradial wave function for (3 — 
{1,1,0,1,1}, J w = 1~ is obtained by multiplying by p 
the wave function for (3 = {0,0,0,0,0}, ,V = 0+. The 
wave function for (3 = {2, 2, 0, 2, 0}, = 2+ is taken to 
be the same as the one for f3 = {2,0,0,0,0}, = 0+. 
It is worth noting that l x + l y < K, while for the ground 



Hamiltonian matrix elements 



Recalling that one set of Jacobi coordinates has been 
selected {x,y,R}, and removing the CM contribution, 
the Hamiltonian is written in hyperspherical coordinates 



H(p,Sl)=T{p,Sl) + V(p,Sl) 



(20) 



The kinetic energy operator, in the differential equations 
for the U(p) = p 5/2 R{p), is 



d 



5 d 



1 



dp 2 pdp p 2 ^ 



(21) 



The kinetic energy operator does not connect channels 
with different (3 or J, i.e. 

(THO; i, (3, J\f\THO; i' , 0', J') = 
(THO; i, (3, J\f\THO; (3, J)Sp,p,Sj,j,, (22) 

and hence only matrix elements within these channels 
have to be calculated. These matrix elements are, with 

4 = s p(p) 2 ' 



j 



(THO; i, P, J\T\THO; i', (3, J) 



N 2 



>'P Lf + 2 (s 2 )L^ 2 (s 2 ) 



dp 



+ ^ dp U o0 (p) 2 Lf_\ 3 (s 2 )L^(s 2 ) (2s {p) 



dsp(p) 
dp 



dUpp(p) \ T K+3( 2^ t K+2( 2n / ^ dsp(p) 

dp p dp 

, dUofi(p) s T K+2, 2\ T K+3, 2n , ^ ds f3 {p) 

o a P a P 



dp ^^L?+ 2 (s 2 )L«+ 2 (s 2 ) 



(23) 



These integrals are evaluated by quadratures. 

The potential energy operator connects different chan- 
nels within the same J. The hyperangular integration is 
performed using the code FaCE 36] on the set of hyper- 



radial quadrature points providing the set of functions 
Vpp, (p). The potential energy matrix elements in the 
THO basis are then simply 



J 



(THO; i,(3, J\V(p, Vt)\THO; i',f3', J) = 



dpU™°{p)V^{p)U™°{p). 



(24) 
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The formalism presented above provides a complete 
basis and the corresponding Hamiltonian matrix ele- 
ments. In the following Section we apply the method 
to study 6 He. 



III. STRUCTURE CALCULATION FOR HE 

The 6 He nucleus is treated here as a three-body sys- 
tem, comprising an inert a core and two valence neu- 
trons. The ground state has total angular momentum 
J n = + and an experimental binding energy of 0.973 
MeV. The ground state wave function was obtained by 
solving the Schrodinger equation in hyperspherical co- 
ordinates, following the procedure described in (2^. l3fil]. 
Besides the two-body (n — n and n — a) potentials, the 
model Hamiltonian also includes a simple central hyper- 
radial three-body force. This is introduced to overcome 
the under-binding caused by the other closed channels, 
the most important of which are t+t channels. The n- 
4 He potential is taken from Ref. H3, with central 
and spin-orbit components, and the GPT NN potential 
with central, spin-orbit and tensor components is 
used. These calculations were performed with the code 
FaCE [36] . All calculations truncate the maximum hy- 
permomentum at K max — 20 and the three-body force 
is adjusted to give the right binding energy. The cal- 
culated three-body wave function has a binding energy 
of 0.954882 MeV and a rms point nucleon matter radius 
of 2.557 fm when assuming an alpha-particle rms matter 
radius of 1.47 fm. 

In Fig.|2]we plot the hyperradial parts for the first three 
channels of the 6 He ground state wave function. The la- 
bels correspond to the quantum numbers {K, l x , l v , I, S}. 
These channels give the main contribution to the wave 
function. 

From the ground state, the LST for the different chan- 
nels are obtained as explained in the preceding section. 
In Fig. Elthe sp{p) for the LST's of the most important 
6 He ground state channels are shown. The THO basis is 
constructed from these LST using Eq. flBt - 

In Fig. 01 the first few hyperradial wave functions for 
the channel [3 — {2, 0, 0, 0, 0}, and for J T = + , are pre- 
sented. This is the most important ground state channel, 
making a 79% contribution to the total norm. We see 
that as the quantum number i increases the wave func- 
tions are more oscillatory and explore larger distances. 

With the THO basis one can proceed to calculate the 
Hamiltonian matrix elements and diagonalize the Hamil- 
tonian matrix in a truncated space. In Fig.03the Hamil- 
tonian eigenvalues for J K = + , for different maximum 
values of the hyperradial excitations, n&, are presented 
up to e = 10 MeV. The calculated ground state energy is 
-0.954886 MeV. 

To study different strength functions from the 6 He 
ground state to the three-body continuum, states with 
different values of J 71 have to be generated. Thus, we 
have to construct a THO basis for these states starting 



— |3=(0,0,0,0,0) 
-- (3=(2,0,0,0,0) 
■-■ |3=(2,1,1,1,1) 




p(fm) 

Figure 2: Radial part for the dominant components of 
the 6 He ground state wavefunction. The labels stand for 
the angular momentum quantum numbers, in the order 

(K,l X ,ly,l,S) 




P(fm) 

Figure 3: Local scale transformations for the first three chan- 
nels included in 6 He ground state wavefunction. 



from the information we have for the ground state. Close 
to the origin the hyperradial wave functions behave as 
p K . Thus, in order to calculate the LST for states with 
,7 = 2 and hypermomentum K (K = 2, 4, 6, ... , 2n) we 
use the ground state component with J = and the same 
value for K (and also l x , l y , I, S if needed). Again, using 
the LST we generate the basis for J = 2 and diagonalize 
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Figure 5: Eigenvalues of J 71 = + states below e = lOMeV. 



the Hamiltonian. 

In Fig. El the eigenvalues of the Hamiltonian for = 
2+ states are presented for different values of rib. Again 
the lowest state is very stable and is close to the energy 
of the known 2 + resonance, that is 0.824 MeV. The states 
with J = 1 have odd K (K = 1, 3, 5, In — 1). Since the 
ground state wave function contains only even K, for the 
generation of the LST for the J = 1 state we have taken 
the ground state component with K — 1 and multiplied 
this by p, to recover the correct behaviour close to the 
origin. To select different channels with the same value 
for K we look for coincidence in l x ,l y ,l,S. Again, from 
the LST we generate the basis for J — 1 and diagonalize 



> 



2 
n, 



4 



Figure 7: Eigenvalues of J = 1 states below e = WMeV . 



the Hamiltonian. 

In Fig. the eigenvalues of the Hamiltonian for J 71 " = 
1~ states are presented for different values of rife. We 
now see that the lowest state changes as a function of rib, 
and goes down in energy as the dimension of the basis is 
increased. 

We now calculate response functions of our system and 
assess their convergence with the dimension of our trun- 
cation of the THO basis. 



A. Completeness 

In general, a pseudostate with energy e n will be a su- 
perposition of the actual continuum states nearby in en- 
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ergy. There are different ways of assig ning an energy 
distribution to a pseudostate [fj Q, E3|- Here we pro- 
pose a method that takes as reference a large basis with 
N t — (Nb + 1) x N c han states, where Nf, is the number 
of hyperradial excitations of the large basis and N c h a n 
is the number of channels for a given J 7r . This basis is 
considered to be complete for the problem under study. 
In this basis an energy distribution /jv(£,£at) is assigned 
to each discrete state (TV = 1, . . . ,Nt). The width of the 
distribution is an increasing function of the energy and 
the distribution can be Gaussian, Lorentzian, or Pois- 
son, etc. We then consider a smaller basis, with nt <C N t 
states, n t = (rib + 1) x N c h a n, where rib is the number of 
hyperradial excitations of the smaller basis. These states 
n (n = 1, . . . , m) can be expanded in the large basis 



J2c(n,N)\N). 



(25) 



TV 



Then, the distribution for the states in the small basis 

/„(£,£„) is 



/„(£,£„) = \(e\n)\ 2 = J2{m^W')C(n,N)C(n,N'r 

NN' 

« J2\(e\N)\ 2 \C(n,N)f 

N 

= ^M£,£ N )|t7(n,iV)| 2 . (26) 



N 



where we have considered that the off-diagonal terms are 
small compared with the diagonal terms in the bigger 
basis. 

We have approximated the /at(£, £jv) by a Poisson dis- 
tribution 



(m + l)( ro+1 ) m ( m + 1 
P{e, en, m) = (rn4 _ u ; e exp | e 



(m+l) i 

e N m\ 
which has a width given by 



t N 



= V(e 2 ) - (£> 2 = 



\Jm, + 1 



(27) 
(28) 
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Figure 8: Completeness of the basis with Nb = 15 and a 
Poisson distribution with m — 20 for the J n = 1~ states. 
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Figure 9: Completeness of the basis with Nt = 15 and a 
Poisson distribution for the J = 2 + states. See text for details 
on the m values used. 



B. Strength functions 



The larger the value selected for m the narrower the dis- 
tribution is. 

In Figs. El and El we present the completeness of the 
large basis to describe states with J 77 = 1~ and J n = 2+ 
respectively. The completeness F(e) is defined as the 
sum over all the excited states of P(e,eN,m). For the 
B{E1) distribution a value of m = 20 was used. For 
the B(E2) the same value of m was used except for the 
region around the resonance for which a value m = 1300 
had to be used in order to get a width consistent with 
the experimental value. It is observed that the basis is 
particularly suited for describing effects in the low lying 
continuum in both cases. 



To study the effectiveness of the THO basis we start 
with global observables related to structure, the strength 
functions. For a given operator O, which couples the 
ground state (n — 1, J 7 * = + ) with excited states with 
angular momentum J w , the following observables are de- 
fined: 

• Total strength 

S T (6,n t ) = ^|(nJ|6|10)| 2 (29) 

71 

where n runs over all the excited states with total 
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Table I: Convergence of different observables as the THO basis 
is increased: calculated ground state energy (eg), expectation 
value in the ground state for p 2 ((p 2 )s), total strength from 
the ground state {St) [this column should converge to the 
column labelled (p 4 )s — (p 2 }% according to the sum rule], the 
ratio between the energy weighted sum and the total strength 
({E(n t ))) and the ratio between the polarizability and the 
total strength ({S _1 (n t ))). 



lib 


e B 


(P 2 ) B 


S T 


( P 4 )b ( P 2 )l 


(E) 


(E- 1 ) 




(MeV) 


(fm 2 ) 


(fm 4 ) 


(fm 4 ) 


(MeV) (MeV -1 ) 





-0.954882 


30.63 


41.6 


855 


27.8 


0.0433 


1 


-0.954884 


30.60 


738 


853 


5.72 


0.194 


2 


-0.954885 


30.62 


851 


856 


3.08 


0.364 


3 


-0.954885 


30.62 


855 


855 


3.01 


0.426 


4 


-0.954886 


30.62 


854 


854 


3.00 


0.433 



Table II: Convergence of observables as a function of 
the dimension (n b ) of the basis considered. The observ- 
ables presented are: the total strength for B(E1) excita- 
tions from the ground state Q2 n -B(-El)io,ni), the energy 
weighted sum (E w ), the ratio between these magnitudes 
((E) = E w /^2 n B(El)io,ni) and the polarizability (a) below 
e = lOMeV. 



n b E„ B ( £1 )i«-«i E w (E) a 





(e 2 fm 2 ) 


(e 2 fm 2 MeV) 


(MeV) 


(e 2 fm 2 MeV 





1.297 


5.38 


4.15 


0.422 


1 


1.323 


5.71 


4.32 


0.403 


2 


1.323 


5.76 


4.36 


0.403 


3 


1.321 


5.76 


4.36 


0.402 


4 


1.320 


5.75 


4.35 


0.402 



angular momentum J w , up to n t = (rib + 1) x N c h a n- 

In the limit of a very large number of states n t , the 
sum rule for the total strength is 

S T (0) = lim S T [d,n t ) 

n t — >oo 

= (io|6 2 |io) - |(io|6|io)| 2 . (30) 

• Energy weighted sum 

nt 

E w (d,n t ) = ^(e„, 7 -e 10 )|(nJ|6|10)| 2 . (31) 

n 

• Polarizability 

a(6,n t ) = 5Z(£„j-eio) _1 |<«J|6|10)| 2 . (32) 

n 

In Table U we present the results for the case O = p 2 
that connects the ground state with the excited states 
with J T = + . This calculation is performed with 
K m ax = 20 and different values of rib. 

Note that, already for rib = 2, the total strength ac- 
quires a value which is very close to the sum rule value. 
This indicates that the THO basis is very efficient to 
describe the ground state through the relatively long 
range monopole operator p 2 . The values of (E(nt)) and 
(_E _1 (n t )) also stabilize for low values of rib and indicate 
the range of energies in the continuum which are relevant. 
In our case, for J 71 = + , it corresponds to low energies 
of around 3 MeV excitation energy above the break-up 
threshold. 



C. B(EA) sum rules 

We will follow the notations and definitions of Brink 
and Satchler The electric multipole operator 



Table III: Same as Table ITT1 but for E2 transitions. 



rib 


E„ -B(-E2)l ,n2 


E w 


(E) 


a 




(e 2 fm 4 ) 


(e 2 fm 4 MeV) 


(MeV) 


(e 2 fm 4 MeV 





3.51 


14.5 


4.09 


1.21 


1 


6.89 


27.6 


4.01 


2.37 


2 


8.41 


31.1 


3.70 


3.06 


3 


8.47 


30.5 


3.60 


3.20 


4 


8.41 


30.1 


3.58 


3.19 



Q\Ms (f) is defined as 



Qxm x (r) = [^jj Z e r x Y XMx (f), (33) 

where e is the electron charge and Z is the atomic number 
of the system. For a nucleus with a core plus two valence 
neutrons, r = y^mfi y /m c is the position of the core 
(charged particle) relative to the center of mass of the 
system. m c is the mass of the core. For 6 He, Z = 2 
and r = y/(2-\/3), as presented in Fig. ^ The reduced 
transition probability is 

B{E\)nj,n' J' = B(EX; nJ -> n'J') 

= \ {n J\\Q x \\n'J')\ 2 (^±^ (34) 

where the reduced matrix element (nJ\ \Q\\ \n' J') is de- 
fined as 

(nJ M\Q\m x \n' J'M') = (-1) 2 A (J'M'AAf A | JM) 

x (nJ\\Q x \\n'J'). (35) 

Using some angular momentum algebra we arrive to the 
final expression 
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OX I 1 

B(E\) nJ>n ,j, = —^\( n J\\Q x \\n'J')\ 2 

47T 



2A 



(2A+l)(2J' + l) z2 g2 / ymJh \ 
47r ' \ m c J 



x ( 'o o' ) W{ - ll ' l y l 'y> Xl ^ W ( JJ ' 11 '' XS ) 



x / / (sin af (cos a) 2 da dp y x U™° (p)* l * l ^ a )^ (a)U™° (p) 



(36) 



Using this expression it is simple to calculate the sum 
rules of electric transitions from the ground state, n = 1, 
,r = 0+, to the states (n, J): Y, n B{E\) 10 , nJ , 

Y,B(E\) w , nJ = r^±lW|(l ||Q A ||nJ)| 2 . (37) 

n ^ ' n 

In particular, we are interested in the A = 1 electric 
dipole strength, connecting the J n = + ground state to 
the J 77 = 1~ states, and in the A = 2 electric quadrupole 
strength, connecting the ground state to the J 71 ' = 2 + 
states: 

• B(E1) sum rule 

2>(M) 10lBl = C% 2 |1 0). (38) 

■^-^ 47T 

• B(E2) sum rule 

^ 5 Z 2 e 2 m 2 u 2 , 

5>(£2) 10 ,„ 2 = j-^ l 0y* 1 0). (39) 

47r rn^ 

i 



where SO(N) = E„(.9 s |0|n)C(n, iV). We have used this 
method to calculate B(EX) distributions. For this pur- 



Tables In] and IIIII show the results obtained for 6 He 
when including states up to 10 MeV in excitation energy 
above the two-neutron break-up threshold. These include 
the B(EX), total strength, energy weighted sum, and po- 
larizability. The B(E1) and B(E2) sum rule values for 
the total strength are 1.500 e 2 fm 2 and 9.99 e 2 fm 4 , respec- 
tively. The values in the second column in Tables ITU and 
IIIII are close to these limits but do not reach them since 
only states up to lOMeV are included. If the complete 
sum is done the values are 1.498 e 2 fm 2 and 9.59 e 2 fm 4 
that are close to the corresponding sum rule values. 



D. B(E1) and B(E2) distributions 

To obtain a continuous B(EX) distribution from the 
discrete values B(E\)io, n j we can apply the procedure 
described in subsection IIII Al In this case we have to 
evaluate 



(40) 



pose we have taken a large basis with N t = (N b + 1) x 
N c han states, where Nb = 15. The Poisson distribution 



|< 5S |6| £ )| 2 = (gs\d\e){e\d\gs) « "£{gs\d\n){n\e){E\n'){n'\d\gs) 

nn' 

= J2(9s\d\n)(n'\d\gs) Y,(K\e)W)C(n,N)C(n',N'r 

nn' NN' 

« J2(9s\d\n)(n'\d\gs)J2C(n,N)f N (e,e N )C(n',Ny 

nn' N 

= ^2fN(e,e N )SO(N)SO(Nr 

N 
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Figure 10: B(E1) distribution for n b — 4 with n t 
1) x N c han and a Poisson distribution with m = 20. 



30 



10 

(n b + 



> 20 



B 10 







_ THOn =4 

b 

- - continuum scattering 



N =15 

b 



2 4 

e (MeV) 

Figure 11: B(E2) distribution for n b — 4 with n t = (n b + 
1) x N chan and a Poisson distribution with m = 20. 



used had m = 20. 

In Figs. QUI and El we present the B(El) and B(E2) 
distributions up to 10 MeV. Calculations for different di- 
mensions of the THO basis are not presented, but are 
practically indistinguishable. For comparison, the dis- 
tribution calculated with the continuum scattering wave 
functions is also shown. The latter is the calculation re- 
ported in Ref. [38]. For the B(E1) the distributions are 
in good agreement with the maximum at around 1.2 MeV 
and with the same total strength. For the B(E2) distri- 
butions the most relevant feature is the appearance of 
a narrow low-lying resonance. In both calculations the 
depth of the three-body interaction was adjusted to re- 
produce the position of the 2 + resonance at 0.824 MeV. 
Again we find good agreement between the calculations 



> 20 
u 



W 10 







THO n =4 full 

b 

THO n =4 resonance 

b 
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_L 
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Figure 12: Total B(E2) distribution for n b = 4 and resonance 
B(E2) distribution for n b — 4 and for a Poisson distribution 
with m = 20. 



in both shape and total strength. As noted above, we 
used a larger m value for the energy region around the 
resonance. Our calculations exhibit a small bump just 
below the resonance which does not appear in the contin- 
uum scattering calculation. This small difference comes 
from the low energy components of the resonance. This 
is more clearly shown in Fig. El where, besides the full 
B(E2) distribution, we have superimposed the contribu- 
tion of the first J 71 ' = 2 + state that appears around the 
resonance energy in the THO calculation for nj, = 4. This 
single state reproduces the resonant peak very well. 



IV. SUMMARY AND CONCLUSIONS 

We have presented a generalization of the transformed 
harmonic oscillator (THO) method proposed in 0] ap- 
propriate for three-body problems. The method provides 
a discrete representation of the continuum spectrum of a 
three-body system from a knowledge of its ground state 
wave function, either in analytical or numerical form. 
This wave function is used to obtain a local scale trans- 
formation that converts this state into a harmonic oscil- 
lator ground state. The inverse transformation is then 
applied to the corresponding excited harmonic oscillator 
states to obtain a basis set for the physical system. Fi- 
nally, the three-body Hamiltonian is diagonalized in this 
basis, providing a discrete representation of eigenstates 
and eigenvalues of the three-body system. 

The formalism has been applied here to the Borromean 
nucleus 6 He, for which several strength functions, in- 
cluding the dipole and quadrupole Coulomb transition 
strengths, have been calculated. These observables are 
found to converge quickly with respect to the number of 
THO basis states included. Furthermore, the calculated 
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strength distributions are in very good agreement with 
previous results obtained using a three-body continuum 
scattering wave functions. The results found in this work 
suggest that the THO basis could also be effective for con- 
tinuum discretization in scattering calculations. Work in 
this direction is in progress and results will be presented 
in a future publication. 
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